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Abstract. We analyze the low-energy e-N2 collisions within the framework of the Modified-Effective Range 
Theory (MERT) for the long-range potentials, developed by O'Malley, Spruch and Rosenberg [Journal of 
Math. Phys. 2, 491 (1961)]. In comparison to the traditional MERT we do not expand the total cross- 
section in the series of the incident momentum hk, but instead we apply the exact analytical solutions of 
the Schrodinger equation for the long-range polarization potential, as proposed in the original formulation 
of O'Malley et at. This extends the applicability of MERT up to few eV regime, as we confirm using 
some simplified model potential of the electron-molecule interaction. The parameters of the effective-range 
expansion (i.e. the scattering length and the effective range) are determined from experimental, integral 
elastic cross sections in the 0.1 - 1.0 eV energy range by fitting procedure. Surprisingly, our treatment 
predicts a shape resonance that appears slightly higher than experimentally well known resonance in the 
total cross section. Agreement with the experimentally observed shape-resonance can be improved by 
assuming the position of the resonance in a given partial wave. Influence of the quadrupole potential on 
resonances is also discussed: we show that it can be disregarded for N2. In conclusion, the modified-effective 
range formalism treating the long-range part of the potential in an exact way, reproduces well both the 
very low-energy behavior of the integral cross section as well as the presence of resonances in the few eV 
range. 

PACS. 34.80.Bm 03.65.Nk 



1 Introduction 

Combining the short and long-term interactions is the 
classical problem of electron-molecule scattering theory. 
The widely used R-matrix method [UEIj for instance, ap- 
plies two distinct methodologies, quantum chemistry and 
long-range interactions for inner and outer regions, respec- 
tively. Other approaches, like that by Czuchaj and collab- 
orators for noble gases [3] use analytic short-range modifi- 
cations to the polarization potential (we skip many other 
similar attempts). Even full ab-initio models use "effective 
core potentials" [4l[5] . Systematic studies of short-range ef- 
fects, including the exact exchange and correlations poten- 
tials have been done for atomic and molecular targets by 
Gianturco and collaborators [HIT]. They showed, both for 
electron and positron scattering how the low-energy 
cross sections are sensitive to a proper choice of the cut-off 
parameters defining the inner and outer regions. 

In our previous paper on positron-atom scattering [lOj 
we have revisited a textbook formulation of the long and 
short range problem: the modified effective range theory 
(MERT) [Tl]. Analyzing the low-energy scattering on Ar 
and N2 targets, we have concluded that MERT expansion 
performed solely on the short-range part of the poten- 
tials, reproduces well the experimental data up to few eV 
regime. Similar approach has been used in [12j . to study 



both electron and positron scattering on CO2, where we 
have shown that MERT, apart from the zero-energy limit, 
is able to reproduce approximately the position of the 
shape resonance at few eV. In the present paper we apply 
MERT to yet another, classical subject of the electron- 
molecule scattering: the N2 molecule in the low-energy 
region, where the well-known shape resonance appears in 
total and vibrational-excitation cross sections (see for ex. 
[l3 l fT4 | fT5]). Shortly, the present approach uses experimen- 
tal elastic cross sections in the very low-energy range (be- 
low 1 eV), see Karwasz et al. [16] to derive the scattering 
potential within the MERT formulation and then calcu- 
lates the integral cross section at energies up to 10 eV. 

In MERT theory the scattering is due to the polar- 
ization potential and the short range potential is mod- 
elled by an effective range expansion. These features make 
MERT particularly suitable for the near-to-zero energy 
range, treated with difhculties in other approaches [TfllTB] . 
The near-to-zero energy MERT expansion continues to be 
applied to numerous problems: integral cross sections for 
atoms [TOll^ . non-polar piU^ and polar molecules 
[24] . MERT is frequently used by experimentalists for ex- 
trapolation of measured differential and integral cross sec- 
tions towards energies and angles inaccessible by experi- 
ment [25 l l26 ] l27] . Different extensions of MERT, consisting 
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in developing phase shifts (or elements of the reactance, 
K matrix) into series of k were applied by Morrison et al. 
[17] to rotational excitation of N2 molecule, and by Macri 
and Barrachina to electron - metastable He scattering and 
electron photodetachment from Li~ 

The most extensive, to our knowledge, analysis of MERT 
applicability in electron scattering was done for noble gases 
by Buckman and Mitroy [30] . They concluded that MERT 
fails above about 1 eV, exact energy depending on the 
target. However, in that as well in numerous other MERT 
applications, those were integral cross sections which were 
developed directly in series of the colliding electron mo- 
mentum k. As we showed in the previous work [lOj . the 
range of applicability of such a series is limited to a few 
tenths of eV. Therein [10] we have investigated an alter- 
native way of applying MERT, that closely follows the 
original formulation of O'Malley et al [TT]. It consists in 
calculating the phase shifts due to the polarization poten- 
tial from Mathieus solutions of the Schrodinger equation 
and introducing the effective-range expansion exlusively 
for the short-range potential. There is no explicit k expan- 
sion for the integral cross sections in our approach; only 
the short-range contribution expressed in terms of an ap- 
propriate boundary condition, is expressed as series of k. 
This makes the numerical procedures more complicated 
but has allowed to extend the applicability of MERT for 
positron-argon and -nitrogen scattering to the range of 2-3 
eV. We note that our method closely resembles quantum- 
defect-thcory formulation of the scattering problem, where 
the quantum-defect parameter, that slowly varies with en- 
ergy, is expanded in a series of k (see e.g. [3T]). In electron- 
atom scattering the quantum-defect theory has been al- 
ready applied by Watanabe and Greene [32] to analyze 
photodetachment process of K~ and by Fabrikant [33] to 
low-energy electron-metal atoms resonances. 

In this paper we apply our MERT-based model to the 
molecule N2 , well studied by beam and swarm techniques 
at energies down to less than 0.1 eV. The parameters of 
MERT expansion for the short-range potential are derived 
from fitting procedure using the recommended [16] inte- 
gral cross sections below 1 cV. Next, we use these poten- 
tials to calculate cross section up to 10 eV, observing a 
broad p-wave resonance around 6 eV. In such an uncon- 
strained fit the resonances appear at higher energies than 
experimental values. By small, within the experimental 
error bar, modifications of the data used for the inver- 
sion procedure, one can reproduce the exact position of 
the resonances. In that case both resonances appearing in 
the p or d-partial wave channels are narrower than exper- 
imentally determined (but this can be partially caused by 
neglecting the nuclear motion, see for ex. Ref. [35] )■ We 
discuss also the relative contributions to integral cross sec- 
tions from the non-spherical part of the polarizability and 
the quadrupole moment of molecule, using the distorted- 
wave approximation. 

The paper is organized as follows. In section [2] we dis- 
cuss the analytical solutions of the Schrodinger equation 
with the polarization potentials, and we present princi- 
ples of the modified effective-range expansion. In section [3] 



we demonstrate that MERT can describe both low-energy 
part and shape resonances in p- and d-wave cross-sections, 
assuming some simplified model of electron-molecule po- 
tential. The actual electron resonance in N2 is analyzed in 
section UJ Section [5] discusses correction to MERT cross 
sections due to the non-spherical part of the polariza- 
tion potential and the molecule quadrupole moment. We 
present some conclusions in section [5] Finally, three ap- 
pendices give some technical details of the analytical so- 
lutions for the polarization potential. 



2 Scattering of a charged particle on a 
molecule and the quantum-defect approach 

The relative motion of a charged particle and a non-polar 
molecule is described by the Schrodinger equation 



2/i 



V{r) - E 



If (r) = 0, 



(1) 



where ^ is the reduced mass, E is the relative energy of 
the particles, and V{y) is the particle-molecule potential. 
Potential V{y) can be written as a sum 



VlT) = Vi{v) + V2{v)+Vs{r) 
of long-range contributions 
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and a short-range potential Vs {r) . The long-range poten- 
tial is expressed in terms of the static spherical and non- 
spherical polarizabilities a and a2, respectively, and the 
quadrupole moment Q. The short-range part Vs(r) comes 
into play at the distances comparable to the size of the 
molecule, where the molecule cannot be treated as a single 
object. Numerous approaches were proposed in the past, 
based on modifications of the polarization potential [5], 
deriving the short-range interaction from electronic densi- 
ties etc. |34j . or expanding the scattering amplitude at low 
energies [33] . In the present approach we express Vs (r) in 
terms of some boundary conditions imposed on the wave 
function at r — > 0, while the effects resulting from the 
finite range of Vs{r) are included explicitly in the frame- 
work of the modified effective-range theory. 

In the following we assume that the nonisotropic part 
V2(r) can be neglected in comparison to the isotropic po- 
larization potential Vi(r). This is justified as long as the 
quadrupole moment Q and the nonsphcrical polarizability 
a2 expressed in atomic units are much smaller than a. In 
such a case the main effect of the non-isotropic terms is 
the coupling between different partial waves, while their 
contribution to the total cross-section at low energies re- 
mains small. As we show later in section Sec. HI where we 
analyze the corrections arising from V2(r), this is a good 
approximation in the case of N2 molecules. 
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The radial part of the Schrodinger equation with the 
isotropic polarization potential Vi (r) reads 



Combining (fTOll . ((34|) . ([35|) one can easily obtain 
the low-energy expansion of the phase shifts [TT] 



_ 1(1 + 1) ^ {R*f ^ 2fi£ 

Qj,2 j,2 J.A f^2 



^i{r)=0, (5) g cot 7/0 



where (r) denotes the radial wave function for the par- 
tial wave /. For the polarization potential, it is conve- 
nient to introduce some characteristic units R* and E*, 
where R* = ^ ae^^/h? denotes the characteristic length 
of the r""* potential, and E* ~ /{2^iR*^) is the charac- 
teristic energy. With appropriate change of variables, the 
Schrodinger equation (O can be transformed into Math- 
ieu's differential equation of the imaginary argument [111 
I351I36| . and solved analytically in terms of the continued 
fractions. Some basic properties of the analytical solutions 
are discussed in Appendices [X| [B| and [C] (see |571[55| for 
more details on Mathicu functions). Here, we focus only 
on the behavior at small and large distances. For r R* , 
when the polarization potential dominates over centrifu- 
gal potential and the constant energy term, behavior of 
^i{r) is given by 



r- , N r^O . I R* , 

Wi (r) ~ r sm h 0/ 

1 J, 
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Here, q = kR* , and a = Us/R* = —I/Aq. The low-energy 
expansions pT|) - P^ are applicable only at small energies 
E ^ E* . In contrast, in our approach we use the exact 
formula ([9]) for the phase shift, with the short-range phase 
expanded according to ([TO]). 



3 Model calculation 

To illustrate how MERT works for the low-energy scatter- 
ing on the long-range potential we consider a simple 
model of the potential with a square-well interaction at 
short distances (see Fig[T]) 



where 4>i is a short-range phase, which is determined by 
Vs(r). For E = and I = the solution ([6]) becomes exact 
at all distances, and from its asymptotic behavior at large 
distances one can easily determine the value of the s-wave 
scattering length 



-R* cot((/)o). 



(7) 



At large distances: r ^ R* , ^i{r) must take the form of 
the scattered wave 



sin(fcr - Z-^ + r]i), 



(8) 



where k = ^/2jlE/h. Using the Mathieu functions one can 
find the following relation between the phase shift rji and 
the short range phase [11] 



tan rji 



— tan(52 -f- tan((/>; + l^)ta.n5{m^ — 1) 
tan (5(1 - m?) + + - tan^ 5) ' 

(9) 

where 5 = ^{v — I — and m and v are some parame- 
ters, that are determined by the analytical solutions (see 
Appendices [^ and |B] for details) . 

In general, the parameter (f>i entering the asymptotic 
formula ([6]), depends on energy, and can be expanded in 
powers of k. In our case it is more convenient to expand 
tan(0i + ), entering formula ([5]): 



tan((/); + ) 



Ai + ^R*Rik^ 



(10) 



where Ai = ta.n{(f>i + )|g_Q- The lowest order correction 
in k is quadratic, and can be interpreted as an effective 
range Ri for the partial wave I [TT] . 



y(r) 




r < R,n 

Rrn <r < Rq 

r>Ro 



(13) 



Here, i?„i denotes a hard-core radius and Rq is the ra- 
dius of the square-well potential. The potential V{r) does 
not reproduce the full behavior of the electron-molecule 
interaction, however, it captures the basic properties of 
shape resonances in such a system, and allows to verify the 
accuracy of our MERT-based approach. The Schrodinger 
equation for the potential V(r) and partial- wave quantum 
number I can be solved exactly: 





f(0 = \ rji{xr)+C{E)ryi{xr) 

A{E)^MJiz) + B{E)^TJ^z) 



where z = In 



r<Rm 
Rm < r < i?o 
r>Ro 

(14) 

E + U, ji{x) 



^/Vkr), ft2^V(2M) = 

and ni{x) denote spherical Besscl functions, and M^{z) 
and T^{z) are the two linearly independent solutions for 
potential. Parameters A{E),B{E), and C{E) are to 
be determined from the continuity conditions for tZ/(r) at 
r = Rni and r = R^. With the help of the asymptotic 
formulas ([30]) and ([3T]) for Mi,{z) and Ty{z) for large and 
negative argument z (large r) we find exact phase shifts 
and the total cross-section for scattering on V{r). 

We compare the exact result with our MERT model, 
where MERT parameters Ai and Ri are determined from 
the expansion of the short range phase Lpi . To this end we 
compare the form of the wave function ([T4]) for r > Rq, 
with the solution ([55]) derived in Appendix B, obtaining 
the following condition for 0;: tan {<j)i{E) + ^v{E) 4- ^) = 
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Fig. 1. The model potential Eq.(IT3} that is used to test MERT 
for the shape resonances 



A{E)/B{E). Expanding oi A{E), B{E) and v{E) in pow- 
ers of k (c.f. (1211)) yields Ai and i?;. 

Fig. 121 shows the dependence of the total elastic cross- 
section on energy, for few particular values of U , for which 
p-wave shape-resonances appear. The range i?o was cho- 
sen to be smaller than i?*, but still, of the same order as 
R* , that roughly corresponds to the conditions of the low- 
energy electron(positron)-molecule scattering. The reso- 
nances show up as peaks at some particular values of the 
energy. Wc observe that MERT, including only two low- 
est expansion coefficients Ai and Ri of the short-range po- 
tential, reproduces very accurately the exact results. The 
accuracy decreases for higher energies, where one can ob- 
serve some discrepancies, in particular small peaks cor- 
responding to the g-wave resonances that appear at dif- 
ferent energies for exact and MERT curves. This results 
from higher sensitivity of the quasi-bound states on the 
parameters of the short-range potential for larger I, and 
the absence of the higher order terms in expansion (|10p 
when E becomes too large. Nevertheless, the approximate 
MERT treatment remains very accurate for the main res- 
onances in p-wave channel, even for energies significantly 
larger than E*. Although, for relatively deep potentials 
considered here, the perfect agreement between the exact 
and MERT results is rather expectable, wc have verified 
that MERT remains quite accurate, even for much weaker 
square well potentials, where the quasi-bound states ex- 
tends well beyond Rq. 

A similar comparison, but for d-wave resonances is pre- 
sented in Fig. l3l In this case the resonance peaks are nar- 
rower due to the stronger centrifugal barrier, and weaker 
coupling between quasi-bound and scattering states. Again, 
the agreement between MERT and the exact results is 
very good, and some discrepancies can be observed only 
at E ^ lO-E*. Finally, we note that for / = 2 and param- 
eters of Fig. 131 the quasi-bound state is totally localized 
within the range of the short-range potential. This shows 
that MERT is able to reproduce well the low energy be- 
havior of the total cross sections, including the regime 
of the shape resonances resulting solely from the short- 
range potential, even if the latter is treated in terms of an 
effective-range expansion. The reason why MERT works 
up to energies much higher than E* can be understood on 
the basis of the shape-independent condition k <^ l/i?o 
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Fig. 2. Total elastic cross-section versus energy for the scat- 
tering in the model potential (|13|l . calculated for Rm = 0.17?*, 
Ro = 0.5R* and for different depths U of the square-well po- 
tential. The exact results for this model potential (solid lines) 
are compared with predictions of MERT (dotted lines, indis- 
tinguishable on the scale of figure) including terms up to for 
the short-range part of the potential. Results are scaled by the 
the characteristic distance R* and the characteristic energy E* 
of the polarization potential. 



1000- 



100 



[7=240.3 



[7=238 



[7=227 




E/E 



Fig. 3. Same as Fig. 12] but for values of U for which shape 
resonances in d wave appear. 



for the short range potential, that estimates the regime 
when the short-range phases 4>i are energy independent. 
For i?o ~ 0.5R*, this leads to the condition E < AE*, 
which is shifted to higher energies, by inclusion of the 
effective-range term ^R*Rik'^ in the expansion (ITOl) . 
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R*{ao) £"(eV) a Ai Ro/R* Ri/R* 

Na 3.397 1.179 0.119 -0.537 0.110 0.262 

Table 1. Characteristic distance R" , characteristic energy 
E* and four fitting parameters: a = a^/R* (s-wave scatter- 
ing length), Ai (zero-energy contribution from the short-range 
potentials for p-wave), Ro (s-wave effective range), and Ri (p- 
wave effective range) in the case of unconstrained fits and using 
only the very low energy experimental data. 



4 MERT analysis of low-energy electron 
scattering 

We investigate the low-energy electron scattering on N2 
molecule using a theoretical model based on MERT for s 
and p partial waves [lOj . Our model contains four unknown 
parameters: the scattering length a = Ug/R* = —I/Aq 
and the effective range Rq for s wave, and the zero-energy 
contribution Ai and the effective range Ri for p wave, 
that are determined by fitting to the experimental data. 
In this way for / ~ and Z = 1 we neglect corrections of the 
order higher than in expansion (|10p . On the other hand, 
for partial waves with / > 2 we include only the lowest- 
order contribution to the scattering phase shifts that are 
due to the long-range polarization potential (fT^ . This is 
sufficient as long as the energy E remains smaller than the 
centrifugal barrier, that for / = 2 has height of 9E*. The 
total clastic cross-section a{k) = 4Trk~^ sin^ r]i{k) is 
determined from the phase shifts ([9]) and , with 1/ and 
m evaluated numerically, using the procedure described in 
Appendices |^ and |B1 

In the calculations we use recent experimental values 
of the polarizability, as measured in electron scattering 
experiments: a = 11.54ao'^ [39]. First, we fit our model by 
the least-square, not weighted method to the experimen- 
tal data exclusively below the resonances, i.e. for energies 
E < 1.2eV, down to the lowest energies available from 
beam experiments |16j . The fitted parameters from this 
check are listed in Table [TJ Note, that the scattering po- 
tential parameters are expressed in characteristic distance 
units. For example, the scattering length amounts to 
0.404 atomic units, close to the value of 0.420 a.u. derived 
by Morrison et al. [17]. Subsequently, we extend MERT 
analysis, using the fitted parameters to higher energies. 
Fig. |4] compare the experimental data (points) with the 
theoretical curves obtained from this fit (solid lines). This 
leads to the resonance maximum in the scattering cross- 
section, that is located at somewhat higher energies that 
the experimental one (at 2.1 eV). This "spontaneously" 
appearing resonance is due to the p-wave, see insert in 
Fig.H 

The range of fitted data corresponds to quite low ener- 
gies in terms of characteristic units: E < E* . Therefore, 
the fitting procedure, within the error bar of experimental 
points, is not able to predict with sufficient accuracy the 
cross section in the energy range corresponding to reso- 
nance peaks. Also, neglection of the nonisotropic part of 
the potential leads to some additional errors. Moreover, 
our scheme neglects the possibility of d-wave or higher 
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Fig. 4. Total cross-section for the scattering of electrons on 
N2 versus the energy. Depicted are: recommended experimen- 
tal data from review [TS] (squares), the theoretical fits based 
on MERT, assuming the resonance (tan 77; = n/2) at 2.1eV in 
s wave (dotted line), p wave (dashed line), and d wave (dot- 
dashed line), and without assumption with respect to the po- 
sition of the resonance (solid line). The inset shows in addition 
the s-wave and p-wave contributions to the MERT fit not as- 
suming the position of resonance. 

I a Al A2 Ro/R* Ri/R* R2/R* 

-0.447 -0.021 — -3.100 -0.104 — 

1 0.157 -4.314 -0.130 -47.05 4.704 0.079 

2 0.123 -0.390 -0.222 -0.301 -0.174 0.293 

Table 2. Six fitting parameters: a = Us/R* (s-wave scattering 
length), Al (zero-energy contribution for p-wave), A2 (zero- 
energy contribution for d-wave), Rq (s-wave effective range), 
-Ri (p-wave effective range) and R2 (d-wave effective range) in 
the case of fits assuming resonances (tan?;; = 7r/2) in s, p and 
d wave. The fits are made to the data below and above the 
experimentally observed resonance (see text for details). Note 
that for s-wave only four parameters (two partial waves) were 
sufficient to fit the data. 



partial wave resonances, that can be suggested e.g. by 
the behavior of differential cross-section for the scattering 
on N2. To eliminate these shortcomings we have subse- 
quently added to the fitting procedure also some experi- 
mental points, located in the energy range above the res- 
onance: 6eV < E < lOeV. In addition we have considered 
three possible scenario assuming that resonance maxima 
occur in partial waves s, p or d, exactly at the energy 
of the resonance peak observed in the experimental data. 
This corresponds to the requirement that the phase shift 
rji assumes value 7r/2 at the peak location for / = 0, 1, 2, 
respectively. In this way we obtain some constrains on the 
six fitting parameters: a, Ai, A2 (zero-energy contribu- 
tions) and Rq, Ri and R2 (effective ranges). The results 
are shown in Fig. lU presenting apart from the previous 
data, also theoretical curves with maxima assumed in s, 
p and d wave. The fitted parameters are listed in Table [2] 
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Wc observe that the height of the resonance peak for s 
wave agrees well with the experimental data (we show ex- 
perimental values averaged over the vibrational structure 
seen in the resonance, following recommended data from 
ref. [IS]). This would suggest a possible s-wave origin of 
the maximum. But this point requires some more detailed 
discussion. 

As said before, the generally accepted classification 
of resonances describes "shape" resonances as occurring 
within the shape of the effective potential, formed from the 
static and polarization contributions with the centrifugal 
barrier added. Obviously, for the s-wave the centrifugal 
barrier does not exist so trapping of the electron within 
the shape of the effective potential is impossible. How- 
ever, our calculations show a quick change in the s-wave 
phase shift, occurring approximately in the region of the 
experimentally observed resonance. On the other hand, 
assumption of resonances in p and d waves would lead to 
much narrower and higher peaks in the total cross section 
than the experimental values. Therefore, a simple inter- 
pretation of the the N2 resonance is not easy: i) elastic 
differential cross sections at the resonance maximum (2.47 
eV) show a dominating p-wave contribution, ii) differential 
cross sections for the first vibrational excitation (zy=0-l) 
at 2.45 eV show almost a perfect d-wave dependence [40] 
while iii) the MERT analysis does not exclude s-wave ori- 
gin of the maximum in the cross-section. Clearly, several 
partial waves contribute to the resonance and our fairly 
simplified MERT analysis is not able to separate these 
contributions. 



5 Corrections due to the nonisotropic part of 
the long-range potential 

Already early works stressed the importance of other, apart 
from the spherical polarization, components of the long- 
range potential in the electron-molecule scattering. For 
N2, non-spherical polarizability is rather strong (see Ta- 
ble[3]), therefore it is necessary to evaluate corrections aris- 
ing from the anisotropic part of the potential. In this sec- 
tion we calculate them, showing that they can be safely 
neglected in comparison to the leading contribution from 
the spherical polarizability. 

In the present analysis we treat the quadrupole and 
non-spherical polarizability contributions together, but us- 
ing two distinct methods for the zero-energy limit and at 
low E ^ E* energies. To analyze the contribution from 
the nonisotropic part of the potential at very low energies 
(i? <C E*), we can compare the small-fc series expansions 
for the whole potential Vi (r) + V2 (r) [51] , to the expansion 
including the spherical part Vi{r) only. For the isotropic 
part, the total cross-section can be obtained from expan- 
sions of the phase shifts given by Eqs. pT]) - P^ . Substi- 
tuting the parameters of the potential listed in the Tabic [3] 
one obtains 



t(ag) Q2(ag) Q(eag) 



N, 



11.54 3.08 



-1.09 



Table 3. Parameters of the long-range potential Vl (r) for Ai'2 
molecule considered in our analysis. For a we use recent exper- 
imental values of |39) . and for 02 and Q we use experimental 
values used in Ref. fTTl. 



On the other hand, application of the small-fc expansion 
for the full potential [Hj leads to 

^IhM. = a^ + 24.93ga + 0.2'Sq + 0.11 

+ 31.63(a - 0.15)(a + 0.16)q^ log(g) + 0{q^) 

(16) 

We observe that for the typical scales of the scattering 
length a ^ 1, the differences in the coefficients of Eq. lfTHl) 
and Eq. p5|) are at most of the order of 10%. 

For higher energies E ^ E* the importance of the 
nonisotropic part can be assessed with the help of the dis- 
torted wave approximation ^42] . This technique is similar 
to the standard Born expansion of the scattering ampli- 
tude, and is based on the assumption that the total po- 
tential can be split into a strong part l^i, that is treated 
exactly, and the remaining contribution V2 , that is treated 
in the pcrturbative manner. In contrast to the Born ap- 
proximation, where the expansion is done around the free- 
space solutions, here the expansion is performed around 
the scattering solution in the presence of Vi . The scatter- 
ing amplitude can be written as a sum of the contributions 
from Vl and V2 separately (see e.g. [43| ) 

/(kf,kO = ((kf|yi|<?+(kO) + (<?-(kf)|T/2|3'+)) , 

(17) 

where (k) is the scattering state for the isotropic part 

^>±(k) = e''*''' + GtVi^^{\^), (18) 
is the free propagator 



= Ymi{E - pV2m ± ie)-\ 



(19) 



gN2(g) 

47r 



24.60ga 31.32g2 \og{q)a'^ + 0{q^) (15) 



and is the scattering state for the full potential, satis- 
fying the following Lippmann-Schwingcr equation 

9+ =$+{'k{) + G+V2^+, (20) 
= lim(£;-pV2m-Fi -(-ie)-\ (21) 

The second term in Eq. (|17p can be written in the form 
of a Born series 

{•^- {U)\V2W+) = (<?-(kf)|F2 + V2GtV2 + . . . |^'+(ki)). 

(22) 

Finally, to calculate the total cross-section we apply the 
optical theorem: a = ^Im/(ki,ki) |42| . 

Taking into account averaging over isotropic distribu- 
tion of the molecule orientations in the sample, one can 
easily show that the lowest order contribution from V2: 
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E/E* gp Ro = 0.5R* Ro = 0.75-R* Rp = R 



1.0 


3.145 


0.0306 


0.0311 


0.0300 


2.0 


4.014 


0.0147 


0.0135 


0.0138 


4.0 


7.428 


0.0059 


0.0052 


0.0053 



Table 4. The total elastic cross-section ctq and corrections 
due to the nonspherical part of the long-range potential in N2 
molecule, for different energies E of the scattered electron and 
different short-range cut-off parameters Ro- 



{<P^ {'k[)\V2\'P^ (ki)) vanishes. Hence, for the small aniso- 
tropic part, the leading contribution is provided by the 
second order term {<1>~ {'kf)\V2Gf V2\'P~^ {ki)) . We evaluate 
this term using analytical solutions for Vi, that both en- 
ter the scattering solution <P='=(k), and the Green function 
Gi, The latter is expressed in terms of partial wave ex- 
pansion including radial solutions for the Vi potential [42] . 
Because of the singular character of the long-range poten- 
tial at r = 0, the second order corrections diverges at 
r — > 0, and we introduce an additional cut-off parameter 
i?o in the integration at short distances. Physically, Rq de- 
scribes the distance when the interaction potential takes a 
form of the long range potential Vi + V2, which is roughly 
of the size of the molecule. 

Table|4]lists the the corrections due to the nonisotropic 
part V2, together with total cross sections coming from 
the isotropic part Vi. The analytical unperturbed solu- 
tion for Vi was calculated for parameters of Table ([T]), 
corresponding to the low-energy unconstrained fits. We 
observe that the corrections are small, of the order of 1%, 
and in principle, they slightly depend on the cut-off pa- 
rameter Rq. The latter is a consequence of the fact that 
the analytical solutions for the pure polarization poten- 
tial Vi become not physical at small distances, where the 
short-range part of the potential comes into play. There- 
fore corrections listed in tables |4] estimate contributions of 
V2 from large distances only, where the potential takes the 
asymptotic form V1-I-V2. On the other hand, we know that 
at short distances the anisotropic contributions in this re- 
gion are not important, since the interaction potential is 
mainly determined by the short-range molecular forces. 



6 Conclusions 

MERT analytical solution was applied previously [10] for 
positron scattering cross sections. In general, positron- 
scattering cross sections in the limit of low energies fall 
monotonically with rising energy, up to the positronium 
formation threshold, and do not exhibit particular struc- 
tures, see [S]. The model [10] approximated well Ar and 
N2 integral elastic cross sections at 0-2 eV. The derived 
scattering potentials were characterized by negative val- 
ues of the scattering length-like parameters: a and Ai, for 
s and p waves, respectively. 

At present, we have applied the MERT analysis to elec- 
tron scattering on N2 molecule. Parameters of the short 
range potential (the scattering length and the effective 
range) for the s, p, and for higher energies also d par- 



tial wave, were obtained via fitting procedures using low- 
energy experimental data. The model approximates well 
the experimental data up to about 1-2 eV and the deduced 
value of the s-wave scattering length agrees well with the 
values known in the literature. 

Deriving the potential parameters from the experimen- 
tal data below resonance and using them again for our 
model at higher energy, we obtained p-wave resonances, 
although at slightly higher energies than observed experi- 
mentally. To check for the consistency of the fit, the model 
was applied jointly to the data below and above resonance, 
assuming that the resonance occur in the s, p or d wave, 
at the energy as observed experimentally. Within the lat- 
ter constraint, a "resonance-like" peak, corresponding to 
the 7r/2 phase change appears also in the s-wave chan- 
nel. In that case the peak is not as sharp as pure shape 
resonances in p and d partial waves. 

As far as amplitude are considered, fixing the reso- 
nance position reproduces well the absolute value of the 
total cross section in its maximum, but only if the reso- 
nance is assumed to appear in the s-wave channel. This 
seems somehow surprising, as the N2 low-energy resonance 
is usually classified as the state. However, the ex- 
tended calculations by Sun et al [13] showed that "elastic 
resonant integral cross sections entail comparable contri- 
butions from the ^Hg and the ^Eg transition matrices [...]" . 
Those authors [l^ showed also that even if a few partial 
waves allow to reproduce the shape of differential cross 
sections, as many as 13 of them are needed in the reso- 
nance to reproduce the absolute values of cross sections. 

More experimental points at low energies, would be 
useful for MERT modelling. We recall here Ramanan and 
Freeman |45j who on the basis of their swarm experiment 
stated: "experimental data presently available do not rule 
out the existence of a Ramsauer-Townsend minimum at 
an electron energy of 0.4 meV or lower" . 

Summarizing, the present work unifies in a single model 
the very low-energy dependence of the integral cross sec- 
tions in N2 with the occurrence of shape resonances at 
a few eV energies. It confirms that limitations on ap- 
plicability of MERT to the very low energy range, come 
rather from direct expansions of the phase-shifts in series 
of fc, rather than from any intrinsic defects of the effec- 
tive range model. Obviously, several aspects of resonant 
scattering, like partitioning between the elastic and vibra- 
tional/rotational excitation are beyond the capacities of 
our model. It would be challenging to incorporate MERT 
potentials to more advanced theoretical treatments, in 
particular of inelastic direct and resonant scattering |461 

mm- 
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A Solutions of the Schrodinger equation for 
potential 

The method for solving the Schrodinger equation with the 
polarization potential was diseussed in our previous paper 
[To] . Below we indicate the basic steps of the derivation. 
In the radial Schrodinger equation (O we substitute r = 
\/^e~^ / ^fk and ^i{r) = ip{r)y/r / R* , which yields the 
Mathieu's modified differential equation [37[|38] 



- [a - 2gcosh2zlV' = 0. 



(23) 



where a = + and q = kR*. Two linearly independent 
solutions M{z) and T{z) can be expressed in the following 
form IMllSZj 

oo 

M,{z) = (-l)"c„(i^)J2„+.(2Vgcoshz), (24) 



n— — oo 



Tu{z) = {-IT Cn{v)Y2n+uM cosh z), (25) 



Here, v denotes the characteristic exponent, and J^(z) 
and Y^(z) are Bessel and Neumann functions respectively. 
Substituting the ansatz ([24]) and ([25|l into ([23|l one obtains 
the following recurrence relation: 

[{2n + vf - a] c„ + 9(c„_i + c„+i) = 0, (26) 

which can be solved in terms of continued fractions. To 
this end we introduce /i+ = c„/c„_i and ft.^ = c_„/c_„+i 
for n > 0, which substituted into (p6|) gives the continued 
fractions 



ht = — 



qh 



+ 

n+l 



d„ ' 



qh, 



(27) 



n+l 



with dn = (2n + v)^ — a. Characteristic exponent has to 
be determined from Eq. ()26|) for n = expressed in terms 
of hi and h^ . In practice, to find numerical values of 
the coefficients c„ we set /i+ = and /i~ = for some, 
sufficiently large m and calculate h^ and h~ up to n = 1 . 
using ((77)) . 



Mu{z) and T^{z), with another pair of solution Afj/(— z) 
and T^{~z) across z — [35]. In this way one obtains ^TU\ 



M,{z) 



Z^ — OQ 



174^ 



TT q 



; (V9e-- + f - f ) (30) 



z/2 



TT q 



1/4 



sm 



cot 7ri/(m^ - 1) cos (y/qe + f - f ) 



(31) 

where m = lim^^Q+ Mi,{z)/M-^{z). To evaluate m nu- 
merically, we rather avoid using formula (|24|) , where sum- 
mation over n converges very slowly for z —f 0. An alter- 
native approach is to use a different representation for the 
solutions of the Mathieu's equation [37] 



W,{z) 



E 

n— — oo 



Cn{v)e 



2n+y 



(32) 



where c„ {v) are the same coefficients as used in definition 
of Mv{z) and Ty{z). One can show that W^{z) differs from 
M,y{z) only by some constant prefactor. We utilize this 
fact calculating the prefactor at some moderate values of 
the argument (z ~ 1), and then evaluate m using the 
function W,y{z), that has well defined behavior at z ^ 0. 

Finally the radial wave function tf'/ (r) for partial wave 
I can be written as linear combination of M^{z) and T^{z) 



Wi{r) =sin(0z + fzy+f) 



R*^ I VR* 
AU In 



cos 




(33) 



where 4>i is a parameter which appear in the small r ex- 
pansion ([6]). Now, the behavior of !f'(r) at small and large 
distance described by Eqs. (H])-®, can be readily obtained 
from asymptotic expansions (|28p - ((3T]) . 



C Expansions for small energies 



B Asymptotic expansions for large arguments 



Behavior of the parameters m and u at small energies is 
given by the following series expansions )49j 



Asymptotic behavior of M^{z) and T^{z) for large z fol- 
lows immediately from asymptotic expansions of Bessel 
functions for large arguments [38j 



M,(z) 



2 e^'^/^ 



(28) 
(29) 



where = Yl^=-oo '^ri{v)- To obtain asymptotic behavior 
for large and negative z it is necessary to connect solutions 



K?) 



and 



m{q) 



+ 0(g4), (34) 



q\i+krC^ 

4 







r(i 







+ 0(g'+^/2), 



(35) 



where ipix) denotes the digamma function [38]. The above 
formulas can be obtained, for instance, by solving the re- 
currence relation (|26p keeping only the lowest-order terms: 
Co, c~i and ci. 
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